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Abstract: Motivated by recent detections of pulsar wind nebulae in very-high-energy (VHE) gamma 
rays, a systematic search for VHE gamma-ray sources associated with energetic pulsars was performed, 
using data obtained with the H.E.S.S. (High Energy Stereoscopic System) instrument. The search for 
VHE gamma-ray sources near the pulsar PSR J1718-3825 revealed the new VHE gamma-ray source 
HESS J1718-385. We report on the results from the HESS data analysis of this source and on possible 
associations with the pulsar and at other wavelengths. We investigate the energy spectrum of HESS 
J1718-385 that shows a clear peak. This is only the second time a VHE gamma-ray spectral maximum 
from a cosmic source was observed, the first being the Vela X pulsar wind nebula. 



Introduction 

It has long been known that pulsars can drive pow- 
erful winds of highly relativistic particles. Con- 
finement of these winds leads to the formation of 
strong shocks, which may accelerate particles to 
^PeV energies. 

The best studied example of a pulsar wind nebula 
(PWN) is the Crab nebula, which exhibits strong 
non-thermal emission across most of the electro- 
magnetic spectrum from radio to >50 TeV 7-rays 
[10]. More recently, VHE 7-ray emission has been 
detected from the VelaX PWN [4], which is an 
order of magnitude older (^11 kyr) than the Crab 
nebula, and its nebula is significantly offset from 
the pulsar position, both in X-rays and VHE 7- 
rays. Offset nebulae in both X-rays and VHE 7- 
rays have also been observed in the Kookaburra 
Complex [2] and for the PWN associated with the 
7-ray source HESS J1825-137 [7, 3]. The latter 
source appears much brighter and more extended 
in VHE 7-rays than in keV X-rays. This suggests 



that searches at TeV energies are a powerful tool 
for detecting PWNe. 

Motivated by these detections, a systematic search 
for VHE 7-ray sources associated with high spin- 
down energy loss rate pulsars was performed, us- 
ing data obtained with the H.E.S.S. instrument. 
The VHE 7-ray data set used in the search includes 
all data used in the H.E.S.S. Galactic plane survey 
[6], an extension of the survey to —60° < I < 
—30°, dedicated observations of Galactic targets 
and re-observations of H.E.S.S. survey sources. It 
spans Galactic longitudes —60° < I < 30° and 
Galactic latitudes —2° < b < 2°, a region cov- 
ered with high sensitivity in the survey. These data 
are being searched for VHE emission from pulsars 
from the Parkes Multibeam Pulsar Survey [11]. 
The search for a possible 7-ray excess is done in 
a circular region with radius 8 = 0.22° (as in 
[6]) around each pulsar position, sufficient to en- 
compass a large fraction of a possible PWN. The 
statistical significance of the resulting associations 
of the VHE 7-ray source with the pulsar is evalu- 
ated by repeating the procedure for randomly gen- 
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erated pulsar samples, modelled after the above- 
mentioned parent population. 

In this search, it is found that pulsars with high 
spin-down energy loss rates are on a statisti- 
cal basis accompanied by VHE emission. The 
search for VHE 7-ray emission near the pul- 
sar PSR J1718-3825 revealed the new VHE 7- 
ray source HESS J 17 1 8 — 385. This paper deals 
with the results from the HESS data analysis of 
HESS J1718— 385 and with its possible associa- 
tions with PSR J 17 1 8 — 3825 and other objects seen 
in radio and X-ray wavelengths. 

H.E.S.S. Observations and Analysis 

The data on HESS J1718— 385 are composed 
primarily from dedicated observations of the 
supernova remnant RX J1713. 7-3946 [1], 
which is located at about 1.6° south-west of 
HESS J1718-385. After passing the H.E.S.S. 
standard data quality criteria based on hard- 
ware and weather conditions, the data set for 
HESS J1718-385 has a total live time of -82 
hours. The standard H.E.S.S. analysis scheme [5] 
is applied to the data, including optical efficiency 
corrections. In this analysis, hard cuts are applied, 
which include a rather tight cut on the shower 
image brightness of 200 photo-electrons and are 
suitable for extended, hard-spectrum sources such 
as PWN. These cuts also improve the angular 
resolution and therefore suppress contamination 
from the nearby RX J1713.7— 3946. To produce 
a sky map, the background at each test position 
in the sky is derived from a ring surrounding this 
position with a mean radius of 1° and a width 
scaled to provide a background area that is about 7 
times larger than the area of the on-source region. 

For spectral studies, only observations in which the 
camera centre is offset by less than 2° from the 
best-fit source position are used to reduce system- 
atic effects due to reconstructed 7-ray directions 
falling close to edge of the field of view. The re- 
maining live time of the data sample is —73 hours. 
The spectral significance is calculated by count- 
ing events within a circle of radius 0.2° from the 
best-fit position, chosen to enclose the whole emis- 
sion region while reducing systematic effects aris- 
ing from morphology assumptions. The proximity 



of the strong source makes it necessary to choose 
the background data from off-source observations 
(matched to the zenith angle and offset distribution 
of the on-source data) instead of from areas in the 
same field of view. For a more detailed description 
of methods for background estimation, see [8]. 

Results 

The detection significance from the search for 
VHE 7-ray emission within 0.22° of the location 
of PSR J1718-3825 is 7.9a. A very conservative 
estimate of the number of trials involved ([6]) leads 
to a corrected significance of 6.2a. 




Right Ascension 

Figure 1 : An image of the VHE 7-ray excess counts of 
HESS J 17 18 — 385, smoothed with a Gaussian of width 
0.06° . The colour scale is set such that the blue/red tran- 
sition occurs at approximately the 3a significance level. 
The black contours are the 4, 5 and 6a significance con- 
tours. The position of the pulsar PSR J1718 — 3825 is 
marked with a green triangle and the Galactic plane is 
shown as a white dotted line. The best-fit position for 
the 7-ray source is marked with a black star and the fit 
ellipse with a dashed line. 

Figure 1 shows the smoothed excess count map 
of the 1° x 1° region around HESS J1718-385. 
A two-dimensional Gaussian brightness profile, 
folded with the H.E.S.S. point-spread function, is 
fit to the distribution before smoothing. Its pa- 
rameters are the width in two dimensions and the 
orientation angle, defined counter-clockwise from 
North. The intrinsic widths (with the effect of 
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the point-spread function removed) for the fit are 
9' ±2' and 4'±1' and the orientation angle is -33°. 
The best-fit position for the centre of the excess is 
RA = 17 h 18 m 7 s ±5 S , Dec = -38°33' ± 2' (epoch 
J2000). H.E.S.S. has a systematic pointing error of 
~ 20". 

For the spectral analysis, a statistical significance 
of 6.8cr (with 343 excess counts) is derived. Figure 
2 shows the measured spectral energy distribution 
for HESS J1718-385 (in E 2 dN/dE representa- 
tion). 



Alternatively, fitting the spectrum by an ex- 
ponentially cut-off power law (dN/dE = 
N E- T e- E ' E ^) gives N Q = (3.0 ± 1.9 stat ± 
0.9 sys ) x lO-^TeV^cm^s- 1 , photon index 
T = 0.7±0.6 sta t±0.2 sys and a cut-off in the spec- 
trum at an energy of E cut = (6 ± 3 sta t ± l sys ) TeV. 
This fit, which is shown as a dashed line in Figure 
2, hasax 2 Mo./. of 1.6/3. 
Both the curved and exponentially cut-off power 
law profiles fit the data well; the former has the 
advantage of showing explicitly the peak energy of 
the spectrum, which has to date only been resolved 
in one other VHE source, VelaX [4]. 
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Figure 2: The energy spectrum of HESS J1718-385, 
which is fit by a curved profile (solid line). Alternatively, 
the fit of an exponentially cut-off power law is shown 
(dashed line, refer to the text for details on both fits). The 
first point in the spectrum lacks statistics due to lower 
exposure at small zenith angles and is plotted as an upper 
limit with at a confidence level of 2a. 

The spectrum is fit by a curved profile (shown as 
the solid line): 
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The peak energy -E pcak is (7 ± l sta t ± l sys ) TeV, 
the differential flux normalisation N = (1.3 ± 
0.3 stat ±0.5 sys ) x 10~ 12 TeV" 1 cm" 2 s" 1 and (3 = 
-0.7 ± 0.3 stat ± 0.4 sys . This fit has a X 2 /d.o.f. 
of 3.2/3. The integral flux between 1 — 10 TeV 
is about 2 % of the flux of the Crab nebula in the 
same energy range [5]. 



Possible Associations 

The 7-ray source HESS J1718— 385 is located 
~0.14° south of the pulsar PSR J1718-3825. 
PSR J1718-3825 appears to be a Vela-like pulsar, 
as it is of comparable age, 90 kyr, and has a similar 
spin period, 75 ms. From the spectral fit of a curved 
profile, the energy flux of HESS J1718-385 be- 
tween (1 - 10) TeV is estimated to 2.9 x 10~ 12 erg 
cm _2 s _1 . With a distance of ^4kpc and a spin- 
down luminosity of E = 1.3 x 10 36 ergs _1 , 
PSR J1718— 3825 is energetic enough to power 
HESS J1718-385, with an implied efficiency of 
e 7 = L 7 /E = 0.5%. 

As can be seen in Figure 3, no obvious X-ray coun- 
terpart is visible for HESS J1718-385. There 
is diffuse extended radio emission, which is par- 
tially coincident with the VHE emission. However, 
this emission seems to be correlated with thermal 
dust emission visible in the IRAS Sky Survey At- 
las [12], suggesting that the radio emission is ther- 
mal and is thus not likely associated with a possi- 
ble PWN. The brightest part of this diffuse feature 
is catalogued as PMN J1717-3846 [14]. From 
the point of view of positional coincidence, ener- 
getics, and lack of other counterparts, the associa- 
tion of HESS J1718 385 with PSR J1718-3825 
seems plausible. To confirm this, additional evi- 
dence from spectral and morphological studies in 
VHE 7-rays and from data at other wavelengths is 
needed. 

HESS J1718-385 may well represent the first 
VHE 7-ray PWN found in a systematic search 
for pulsar associations, despite the present lack 
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HESS J1718-385 
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Figure 3: Radio image from the Molonglo Galactic 
Plane Survey at 843 MHz [9] (in Jy/beam). The H.E.S.S. 
significance contours are overlaid in black and the pul- 
sar position is marked with a red triangle. Adaptively 
smoothed ROSAT hard-band X-ray contours are shown 
in green [13]. 

of a PWN detection in other wave bands. The 
remarkable similarity between HESS J1718— 385 
and other known VHE PWNe, together with the 
lack of other probable counterparts, gives addi- 
tional confidence. The detection of an X-ray PWN 
would provide confirmation. 
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